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Abstract. Quantum dynamics (i.e., the Schrodingcr equation) and classical 
dynamics (i.e., Hamilton equations) can both be formulated in equal geometric terms: 
a Poisson bracket defined on a manifold. In this paper we first show that the 
hybrid quantum-classical dynamics prescribed by the Ehrenfest equations can also 
be formulated within this general framework, what has been used in the literature 
to construct propagation schemes for Ehrenfest dynamics. Then, the existence of a 
well defined Poisson bracket allows to arrive to a Liouville equation for a statistical 
ensemble of Ehrenfest systems. The study of a generic toy model shows that the 
evolution produced by Ehrenfest dynamics is ergodic and therefore the only constants 
of motion are functions of the Hamiltonian. The emergence of the canonical ensemble 
characterized by the Boltzmann distribution follows after an appropriate application 
of the principle of equal a priori probabilities to this case. Once we know the canonical 
distribution of a Ehrenfest system, it is straightforward to extend the formalism of 
Nose (invented to do constant temperature Molecular Dynamics by a non-stochastic 
method) to our Ehrenfest formalism. This work also provides the basis for extending 
stochastic methods to Ehrenfest dynamics. 
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1. Introduction 

The Schrodinger equation for a combined system of electrons and nuclei enables us 
to predict most of the chemistry and molecular physics that surrounds us, including 
bio-physical processes of great complexity. Unfortunately, this task is not possible in 
general, and approximations need to be made; one of the most important and successful 
being the classical approximation for a number of the particles. Mixed quantum-classical 
dynamical (MQCD) models are therefore necessary and widely used. 

The so-called "Ehrenfest equations" (EE) result from a straightforward application 
of the classical limit to a portion of the particles of a full quantum system, and constitute 
the most evident MQCD model, as well as a first step in the intricated problem of mixing 
quantum and classical dynamics. For example, it can be noted that much of the field of 
molecular dynamics (MD) -whether ab initio or not- is based on equations that can be 
obtained from further approximations to the Ehrenfest model: for instance, written in 
the adiabatic basis, the Ehrenfest dynamics collapse into Born-Oppenheimer MD if we 
assume the non-adiabatic couplings to be negligible. 

It is not the purpose of this work, however, to dwell into the justification or validity 
conditions of the EE (see for example Refs. [H El El SI [5] for rigurous analyses). Nor is 
it to investigate the unsettled problem of which is the best manner of mixing quantum 
and classical degrees of freedom. The Ehrenfest model has a proven application niche, 
and for this reason we are interested in investigating some of its theoretical foundations, 
in a manner and with an aim that we describe in the following. For recent progress in 
non-adiabatic electronic dynamics in MQCD see, for example, Ref. [6] 

Classical mechanics (CM) can be formulated in several mathematical frameworks, 
each corresponding to a different level of abstraction (Newton equations, the 
Hamiltonian formalism, the Poisson brackets, etc.). Perhaps its more abstract and 
general formulation is geometrical, in terms of Poisson manifolds. Similarly, quantum 
mechanics (QM) can be formulated in different ways, some of which resemble its classical 
counterpart. For example, the observables (self-adjoint linear operators) are endowed 
with a Poisson algebra almost equal to the one that characterizes the dynamical variables 
in CM. Moreover, Schrodinger equation can be recast into Hamiltonian equations 
form (see [7]) by transforming the complex Hilbert space into a real one of double 
dimension; the observables are also transformed into dynamical functions in this new 
phase space, in analogy with the classical case. Finally, a Poisson bracket formulation 
has also been established for QM, which permits to classify both the classical and the 
quantum dynamics under the same heading. 

This variety of formulations does not emerge from academic caprice; the successive 
abstractions simplify further developments of the theory, such as the step from 
microscopic dynamics to statistical dynamics: the derivation of Liouville equation (or 
von Neumann equation in the quantum case), at the heart of statistical dynamics, is 
based on the properties of the Poisson algebra. 

The issue regarding what is the correct equilibrium distribution of a mixed 
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quantum-classical system is seen as a very relevant one [El |9], HDl ttH [iSl [13], [HI [IHl [IE] . 
An attempt to its derivation can be found in Ref. [H], where they arrive to the same 
distribution that we will advocate below, although it is found using the Nose-Hoover 
technique pTl [18] . which in principle is only a mathematical scheme to produce the 
equilibrium, and which the very authors of Ref. [2] agree that it is not clear how to 
apply to a mixed quantum-classical dynamics. On the other hand, in Ref. [8] they 
provide some analytical results about this distribution, but not in the case in which the 
system of interest is described by the EE; instead (and as in Refs. [91110]), they assume 
that the system is fully quantum, and that it is coupled to an infinite classical bath via 
an Ehrenfest-like interaction. 

It is therefore necessary to base Ehrenfest dynamics -or any other MQCD model- 
on firm theoretical ground. In particular, we are interested in establishing a clear path 
to statistical mechanics for Ehrenfest systems, which in our opinion should be done by 
first embedding this dynamics into the same theoretical framework used in the pure 
classical or quantum cases (i.e., Poisson brackets, symplectic forms, etc). Then the 
study of their statistics will follow the usual steps for purely classical or purely quantum 
ensembles. In this respect, it should be noted that other approaches to MQCD (not 
based on Ehrenfest equations) exist [HI [20l [21], and to their corresponding statistics 
[221 123] , however it has been found that the formulation of well defined quantum-classical 
brackets (i.e., satisfying the Jacobi identity and Leibniz derivation rule) is a difficult 
issue [211 [23 [2ni I2Z] ■ On the contrary, as shown below, and as a result of the fact that 
the evolution of a quantum system can be formulated in terms of Hamilton equations 
similar to those of a classical system [H [21] , we will not have difficulties in a rigorous 
formulation of the Ehrenfest dynamics in terms of Poisson brackets. 

The roadmap of this project is the following: In Section [2] we recall the definition 
of the Ehrenfest model. In Section [3] we quickly summarize the formulation of CM 
in terms of Poisson brackets. Then we summarize, in Section |4| the analogous 
description of QM in terms of geometrical objects which can be found, for example, 
in Refs. [2Sl[71[3ni[2Il[321[331l31l3Sl[3S]. These works demonstrate how Schrodinger 
equation can be written as a set of (apparently classical) Hamiltonian equations. Also, 
by using a suitable definition of observables as functions on the (real) set of physical 
states, Schrodinger equation can be written in terms of a canonical Poisson bracket. QM 
appears in this way as "classical", although the existence of extra algebraic structure 
encodes the probabilistic interpretation of measurements, and the superposition and 
indetermination principles [23 [301 I3ll [321 |33]. As we have in mind the application 
to computer simulations, we will consider only finite dimensional quantum systems, 
obtained by a suitable sample of the electronic states. This simplifies significantly the 
description, even if the framework is valid also for the infinite dimensional case (see 
[33]). 

For the case of Ehrenfest dynamics, we can now combine it (Section |5]) with 
CM, whose Poisson bracket formulation is very well known. This is the first of the 
contributions of the paper: once we realize that the geometric description of classical 
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and quantum dynamics are formally analogous, we can combine them as we do when 
combining different classical systems. The procedure is very simple: we consider as 
a global phase space the Cartesian product of the electronic and the nuclear phase 
spaces and define a global Poisson bracket simply as the sum of the classical and the 
quantum ones. It is completely straightforward to prove that such a Poisson bracket is 
well defined and that it provides the correct dynamical equations (EE). From a formal 
point of view, the resulting dynamics is more similar to a classical system than to a 
quantum one, although when considering a pure quantum system the dynamics is the 
usual Schrodinger one. This is not surprising since the coupling of the classical and 
quantum systems makes the total system nonlinear in its evolution, and this is one 
of the most remarkable differences between a classical and a quantum system in any 
formulation. 

Once the dynamical description as a Poisson system is at our disposal, it is a very 
simple task (Section [6]) to construct the corresponding statistical description following 
the lines of Refs. [37] and [38]. This is another main result of the paper: not having 
proved the Hamiltonian nature of Ehrenfest dynamics, it was not possible to provide 
a self-consistent definition of the Statistics corresponding to it. Once the Hamiltonian 
description is available, the definition of this Statistics is straightforward. 

The formal similarity with the CM case ensures the correctness of the procedure 
and allows us to derive a Liouville equation. The choice of the equilibrium distribution 
is based, as usual, in the principle of equal a priori probablities implicitely used by 
Gibbs and clearly formulated by Tolman [39]. After this, in Section [t] we extend Nose 
formalism to the Ehrenfest MQCD framework and provide then a method to simulate 
the canonical ensemble of Ehrenfest systems. This is the final main contribution of the 
paper. 

Please note that for systems where the Born-Oppenheimer approximation is not 
sufficient, the use of the canonical ensemble of a mixed quantum-classical system such 
as Ehrenfest 's is compulsory, since those systems where the electronic gap is small can 
be affected by electronic temperature effects of the same order. Therefore, a simple 
mechanism of simulation such as Nose is very important for practical applications. 

Finally, in the conclusions (Section [s]), we propose to extend, within our formalism, 
stochastic methods to Ehrenfest dynamics. 

2. The Ehrenfest model 

The Ehrenfest Equations have the following general form: 




(1) 



(2) 



(3) 
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where [R, P) denote collectively the set of canonical position and momenta coordinates 
of a set of classical particles, whereas ip is the wavefunction of the quantum part of the 
system (see Ref |10] for the issue of the Hellmann-Feynman theorem in this context, 
i.e., whether one should take the derivative inside or outside the expectation value). 

This and other MQCD models appear in different contexts; in many situations, 
the division into quantum and classical particles is made after the electrons have been 
integrated out, and is used to quantize a few of the nuclear degrees of freedom. However, 
a very obvious case to use EE is when we want to treat electrons quantum mechanically, 
and nuclei classically. In order to fix ideas, let us use this case as an example: the 
Hamiltonian for the full quantum system is: 

1 o , o \ ^ 1 o 1 \ ^ Z jZx 



H :-- 



1^1 1 



4vreo^|r-; -f,| 4neo \Rj - f,\ 



2M 

(4) 

J 

where all sums must be understood as running over the whole natural set for each 
index. Mj is the mass of the J-th nucleus in units of the electron mass, and Zj is the 
charge of the J-th nucleus in units of (minus) the electron charge. Also note that we 
have defined the nuclei-electrons potential Vn-eir,R) and the electronic Hamiltonian 
He{f,R) operators. 

The EE may then be reached in the following way [H El |3l HI [5]: first, the full 
wave function is split into a product of nuclear and electronic wave functions, which 
leads to the time dependent self-consistent field model, in which the two subsystems are 
quantum and coupled. Afterwards, a classical limit procedure is applied to the nuclear 
subsystem, and the EE emerge naturally. In terms of the nuclei positions Rj and of the 
element of the Hilbert space {ip) G "H which encodes the state of the electrons of the 
system, the Ehrenfest dynamics is then given by 

MjRj =-{^\VjHe{f,Rm, (5) 

thj^\i;) = H,ir,RM. (6) 

These equations can be given a Hamiltonian-type description by introducing a 
Hamiltonian function of the form: 

H{R,p) = E ^ + {me{f,Rm- (7) 



Then, by fixing a relation of the form Pj = MRj, we obtain a structure similar to 
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Hamilton equations: 



= (8) 



J 



Pj =-{ilj\WjH,{r,Rm, (9) 

thj^\i;) = H,{r,Rm. (10) 

This set of equations perhaps constitutes the most evident MQCD model, as well 
as a first step in the intricated problem of mixing quantum and classical dynamics. In 
this sense, it is interesting to note that written in the adiabatic basis {ipm{.R)}m formed 
by the eigenvectors of the electronic Hamiltonian, 

He{f,RMm{R)) = E^RMmiR)) , (n) 

the Ehrenfest dynamics collapse into ground state Born-Oppenheimer MD (gsBOMD) if 
we assume the non-adiabatic couplings to be negligible. Indeed, by expanding the wave 
function in this basis (some relevant considerations related to this change of variables 
can be found in Ref . |10] , Section 2 "Ehrenfest Dynamics: Fundaments and implications 
for First Principles Simulations"), 

m 

the equations are transformed into: 

MjRjit) = -J2 \cmit)\^^jEURit)) 

m 

- E cit)cn{t) [Eum) - EMt))] drm)) . (is) 



d ^ 

i-T:Cm{t) = Em{R{t))Cm(t) - 2 ^ C„(t) 



(14) 



dt 

n 

where the "non-adiabatic couplings" are defined as: 

dT{R) ■■= {^m{R)\^jMR)) ■ (15) 

If these are negligible, and we assume that the electronic system starts from the ground 
state (cm(0) = Smo), the EE model reduces to gsBOMD: 

Mj'Rj{t) = VjEom)) , (16) 

Cm{t) =SmO- (17) 

In spite of the formal similarities, the EE equations do not correspond yet to 
Hamilton equations, since they lack a global phase space formulation (encompassing 
both the nuclear and the electronic degrees of freedom) and a Poisson bracket. 
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3. Classical mechanics in terms of Poisson brackets 

We begin by recalling very quickly the Hamiltonian formulation of classical dynamics 
(we address the interested reader to a classical text, such as Ref. [H], for a detailed 
presentation). Let us consider a classical system with phase space Mc, which, for the 
sake of simplicity, we can identify with M^", where n is the number of degrees of freedom 
(strictly speaking, Mc is a general 2n-dimensional manifold, homeomorphic to M^" only 
locally). The 2n dimensions correspond with the n position coordinates that specify 
the configuration of the system, and the n corresponding momenta (mathematically, 
however, the division into "position" and "momenta" coordinates is a consequence of 
Darboux theorem - see below). 

The "observables" in classical mechanics are differentiable functions 

f:Mc^R, (18) 

that assign the result of a measurement to every point in Mc- On this set of functions 
C°°{Mc) we introduce the Poisson bracket, {-, ■}, a bilinear operation 

{-, ■} : C^{Mc) X C^{Mc) ^ C^{Mc), (19) 

which: 

• is antisymmetric, 

{f,9} = -{9j}, yf,geC^{Mc), 

• satisfies the Jacobi identity, i.e., V/, h G C°°{Mc): 

{/, {9,h}} + {h,{f,g}} + {g,{hJ}} = 0, 

• and satisfies the Leibniz rule ,i.e., Wf,g, h G C°°(Mc): 

{f,gh} = {f,g}h + g{f,h}. 

If the Poisson bracket is non degenerate (it has no Casimir functions), a theorem due 
to Darboux ensures that there exists a set of coordinates {R, P) for which the bracket 
has the "standard" form (at least locally, in a neighbourhood of every point): 

If fl = V^^_^^ (on) 

t/l'/2| l^Qp^Q^k QRkQp^- ^^^^ 

k=l 

These coordinates are called Darboux coordinates. They are specially useful when 
studying the dynamics and invariant measures of a Hamiltonian system [H] . In the rest 
of the paper we shall be working with this kind of coordinates. 

The Poisson bracket allows us to introduce the concept of Hamiltonian vector field: 
Given a function / G C^{Mc) and a Poisson bracket {■, ■}, a vector field, Xj, is said to 
be its Hamiltonian vector field if 



Xf{g) = {f,g}, \fgEC^{M, 



c) 
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In standard coordinates: 

^ _^ df{R,P) d df{R,P) d 
^ ^ dPk dW" dR^ dPk 

k=l 

We shall call a Hamiltonian system to a triple (Mc {-,■}, iJ), where {■, ■} is 
a Poisson bracket on Mc, and H G C°°(Mc) is the Hamiltonian of the sytem. The 
dynamics of a Hamiltonian system can be formulated in two alternative manners: 

• The trajectories of the system are given by the integral curves of the Hamiltonian 
vector field Xh- 

• If we consider instead the set of classical observables, the dynamics is written as 
the Poisson bracket of the Hamiltonian function H with any other function, i.e., 

| = {if,/}, V/gC°°(Mc). (22) 

Both approaches are equivalent: the differential equations that determine the integral 



curves of the Hamiltonian vector field are given by Eq. (22), for the functions "position" 
and "momenta" of each particle (i.e., R and P). These equations are nothing else but 
Hamilton equations: 



p. = - ^ (24) 



4. Summary of geometric quantum mechanics 

The aim of this section is to provide a description of quantum mechanical systems by 
using the same geometric tools which are used to describe classical mechanical systems, 
outlined above. It is just a very quick summary of the framework which has been 
developed in the last 30 years and which can be found in Refs. [2S1 El EOl EH E21 ESI 
EH ESI ES] and references therein. For the sake of simplicity, we shall focus only on the 
finite dimensional case. The Hilbert space "H becomes then isomorphic to for n a 
natural number, ^ 7i: 

4.1. The states 

Consider a basis {|V'a:)} for H. Each state {ip) G can be written in that basis with 
complex components {zk}: 

k 

We can just take the vector space inherent to the Hilbert space, and turn it into a real 
vector space Mq, by splitting each coordinate in its real and imaginary part: 

Zk = qk + iPk ^ iqk,Pk) e M^" = Mq. 
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We will use real coordinates {qk,Pk), k = 1, . . . ,n to represent the points of H when 
thought as real manifold elements. From this real point of view the similarities between 
the quantum dynamics and the classical one described in the previous section will be 
more evident. Sometimes it will be useful to maintain the complex notation or Zk for 
the elements of the Hilbert space. Please notice that, despite the use of "q" and "p", 
these coordinates have in principle no relation at all, with the position or the momentum 
of the quantum system. They are simply the real and imaginary parts of the complex 
coordinates used to represent the (finite-dimensional) Hilbert space vector in the chosen 
basis. Therefore, the resulting geometrical description has no relation with the usual 
phase-space descriptions of Quantum Mechanics, such as the Wigner representation, 
which is defined for infinite-dimensional quantum systems and where the coordinates to 
represent the point of the phase space do correspond to the spectra of the position and 
momentum operators. For the same reason, there is no relation with any classical limit 
of the quantum system: the description in terms of these coordinates {q,p) is purely 
quantum. 

There is no one-to-one correspondence between physical states and elements of 
the Hilbert space used to describe it; it is a well known fact that physical states are 
independent of global phases. This ambiguity also translates to the formulation with 
real vector spaces. Two equivalent states are related by a phase transformation; the best 
way to characterize these transformations in Mq is by introducing their infinitesimal 
generator, given by: 



k 



This is a vector field whose meaning is clear if we realize that a phase change modifies the 
angle of the complex number representing the state, when considered in polar form (i.e.. 



in polar coordinates {rk,Ok}k=i,---,n with = r^e Eq. (25) becomes F = 



Then, from a geometrical point of view we can use Eq. (25) in two ways: 

• Computing its integral curves, which are the different states which are obtained 
from an initial one by a global phase multiplication. 

• Acting with the vector field on functions of Mq. A function will be acceptable as 
an observable if it is invariant under phase transformations, i.e., F/ = 0. 

In a quantum Hilbert space, we also consider two states to be equivalent if they 
merely differ by their norm. One may therefore only consider the sphere of states with 
norm equal to one in C"; this corresponds in the real-vector-space description to the 
{2Nq — 1) -dimensional sphere: 



SQ:=\{q,p)eMQ\Y,iql+Pl 

I k 



i+K) = lj- (26) 

It is immediate that the vector field F is tangent to Sq, since the phase change preserves 
the norm of the state. 
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4-2. The Poisson bracket 

Taking as coordinates = qk + ip^, the components of a vector in an orthonormal basis 
we can define a Poisson bracket in Mq, by: 



4^ 2 \dpkdqk dqkdpk) 



(27) 

k 

that corresponds to the standard Poisson bracket in classical mechanics. It will also be 
useful to introduce a symmetric bracket by 

df dg , df dg 



dqk dpk dpk 



(2J 



We have made use of a specific basis; one may think that the definition of the brackets 
could depend on the choice of the basis. However, it is a matter of simple algebra to 
prove that this is not the case. 

An important property of the brackets defined above is that they are preserved by 



the vector field P in Eq. (25) in the sense that 

T{f,g} = {Tf,g} + {f,Tg} 

and 

r{f,gU = {rf,g}^ + {f,rgU. 

A fact that can be proved with a simple computation. This property is important for us 
because it implies that the symmetric or antisymmetric bracket of two phase invariant 
functions will also be phase invariant: 

Tf = Tg = O^T{f,g} = T{f,g}+ = 0. 



4-3. The observables 

We now proceed to discuss how to represent the physical observables on this new setting. 
Instead of considering the observables as linear operators (plus the usual requirements, 
self-adjointness, boundedness, etc.) on the Hilbert space "H, we shall be representing 
them as functions defined on the real manifold Mq, as it is done in classical mechanics. 
But we can not forget the linearity of the operators, and thus not any function is 
acceptable. We will consider, for any operator A G Lin('H), the quadratic function 

fA{^) := imi^)- (29) 

We shall denote the set of such functions as J-'s(Mq). 

Notice that this definition of observable is different from the analogous one in the 
classical case. In the classical framework a state, represented by a point in Mc, provides 
a well defined result for any observable / : Mc — M. In this geometric quantum 
mechanics (GQM in the following) framework, on the other hand, the value at a given 
state of fA G J^s{Mq) provides just the average value of the corresponding operator in 
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that state. Besides, contrarily to the classical case, not all the C°° functions on Mq are 
regarded as observables, but only those of the form of (29). 

Once the definition has been stated, one must verify that it is a consistent one, 
meaning that (i) they are phase invariant functions (F/ = 0), and (ii) the set of 
observables is closed with respect to both the symmetric and antisymmetric brackets 
(i.e. the bracket of two elements of J-'s(Mq) is also an element of J-'s{Mq)). Both 
facts can be proven by direct computation. The second proof, however, is interesting 
because it leads us to consider the relevant algebraic operations within J^j(Mq). In the 
usual approach of quantum mechanics, there are three algebraic structures on the set of 
operators, which turn out to be meaningful and important for the physical description: 

• The associative product of two operators: 

A,B e Lin(H) ^A-Be Lin(H). 

It is important to notice, though, that this operation is not internal in the set 
of Hermitian operators (i.e., those associated to physical magnitudes), since the 
product of two Hermitian operators is not Hermitian, in general. 

• The anticommutator of two operators: 

A, Be Un{H) [A,B]+ G Un{n). 

• The commutator of two operators: 

A,B e Lin(H) ^ i[A,B] G Lin(H). 

Notice that these last two operations are internal in the space of Hermitian operators. 
How do we translate these operations into the GQM scheme? If we take fAifs ^ 

-^s(Mq), 

• the anticommutator of two operators, becomes the symmetric bracket or Jordan 
product of the functions (see Ref. [42]): 

• The commutator of two operators, transkates into the Poisson bracket of the 
functions: 

i[A,B] ^ f\yAM = {fAjs} G J-.(Mq). 

We may conclude, therefore, that the operators and their algebraic structures are 
encoded in the set of functions with the brackets associated to the Hermitian product. 
Thus we see why it makes sense to consider only that specific type of functions: it is a 
choice which guarantees to maintain all the algebraic structures which are required in 
the quantum description. 

Another important property of the set of operators of quantum mechanics is the 
corresponding spectral theory. In any quantum system, it is of the utmost importance 
to be able to find eigenvalues and eigenvectors. In GQM, these objects are considered 
in the following way: Let be the function associated to the observable A. Then, if 
we consider ip E Sq, 
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• the eigenvectors of the operator A coincide with the critical points of the function 
fA, i.e., dfAi'ip) = =r' '0 is an eigenvector of A. 

• The eigenvalue of A at the eigenvector ip is the value that the function takes at 
the critical point ip. 

We finalize this discussion on the observables by noting that the identity operator 
I corresponds with the function 

Mq,p)-J2(^l+Pl)^ (30) 

k 

and the Hamiltonian vector field of this function is precisely the generator of phase 
transformations : 

{/..■}-E(-|--*|;)-r. (31) 

4.4- The dynamics 

As in the classical case, the dynamics can be implemented in different forms, always in 
a way which is compatible with the geometric structures introduced so far: 

• In the Schrodinger picture, the dynamics is described as the integral curves of the 
vector field Xf^, where H is the Hamiltonian operator of the system: 

Xf, = h-'{fH, •}. (32) 

• In the Heisenberg picture, the dynamics is introduced by transferring the Heisenberg 
equation into the language of functions: 

fA^h-'{fH,fA}. (33) 

For example, the von Neumann equation is rewritten as: 

fp = fi-'{fHjp}, (34) 

where fp is the function associated to the density matrix. 

Therefore, it is possible to describe quantum dynamics as a flow of a vector field on the 
manifold Mq or on the set of quadratic functions of the manifold. Such a vector field is 
a Hamiltonian vector field with respect to the Poisson bracket encoded in the Hermitian 
structure of the system. 

The dynamics thus defined preserves the norm of the state, which allows us to 
restrict the set of admissible states to those with norm equal to one. This can be proven 
by analyzing the evolution of the function fi, given by the Poisson bracket with the 
Hamiltonian: 

fi = h-'{fHji}. (35) 

But because of the properties of the bracket, {fn, fi} = f[H,i\ = 0, and thus the flow is 
restricted to the sphere Sq. 
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Example 1 Let us consider a simple two-level quantum system, defined on C^. As a 
real manifold, Mn ~ M^. Consider then a Hamiltonian if : — ?■ C^; 



H 



H21 



H12 
H22 



(36) 



// we consider it as a matrix on the real vector space Mq, it reads: 







-^51 92 


Hqip2 








Hp^P2 


-^92 91 




-^9292 


^q2P2 


Hp2qi 


Hp2pi 


Hp2q2 


Hp2P2 



(37) 



It is easy to prove that this matrix is symmetric if the Hamiltonian is Hermitian 
(H12 = H*^). 

The function fn in J^s{Mq) becomes thus: 

fn = V'K-fiK^iR, where = (gi,Pi, '?2,P2)* 
and then, the Hamiltonian vector field turns out to be: 

'dfn d dfn d , dfn d dfn d 



X 



H 



h- 



dpi dqi dqi dpi dp2 dq2 dq2 dp2 



and its integral curves are precisely the expression of Schrddinger equation when we 
write it back in complex terms: 

^^^{Hp^q^qi + Hp^q2q2 + Hp^p^pi + Hp 



^PlP2 



P2) 



Pi = {Hq^q^qi + Hq^q2q2 + Hq^p^Pi + Hg-^p2P2), 

q2 = h ^{Hp2q2q2 + Hp2q^qi + Hp2p2P2 + Hp2p^pi), 

P2 = -^''^{Hg2q2q2 + Hq2q^qi + Hg2p2P2 + -f^qapiPl)- 

We can write these equations as: 



with 




1 





-1 










1 








(3J 



or, equivalently, 



-h-HH\ 



where 



qi + ipi 
q2 + ip2 



(39) 



is the complex state vector in terms of the real coordinates. The operator J, that satisfies 
= —I , is called the complex structure. As it is apparent when comparing (38) and 
(39), it implements the multiplication by the imaginary unity i in the real presentation 
of the Hilbert space. Observe that J and H^ commute, this is due to the fact that the 
latter comes from an operator H in the complex Hilbert space. 
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5. Ehrenfest dynamics as a Hamiltonian system 

In this section we show how to put together the dynamics of a quantum and a classical 
system following the presentation of the two previous sections. We are describing thus 
a physical system characterized by the following elements: 



5.1. The set of states of our system 

• First, let H be a Hilbert space which describes the quantum degrees of freedom of 
our system. For example, it could describe the electronic subsystem; in this case, it 
is the vector space corresponding to the completely antisymmetric representation 
of the permutation group Sjy (i.e., a set of Slater determinants), where N is the 
number of electrons of the system and each electron lives in a Hilbert space of 
dimension M. Thus, the dimension of H will be Nq = j^T(J[z^,- 
We know that it is a complex vector space, but we prefer to consider it as a real 
vector space with the double of degrees of freedom and denote it as Mq. Also, in 
correspondence with the Hilbert space vectors in the usual formalism of quantum 
mechanics, several states in Mq represent the same physical state. To consider true 
physical states one should extract only those corresponding to the projective space, 
which can be identified with a submanifold of Mq. A more general approach is to 
consider the sphere of states with norm equal to one, Sq, and take into account the 



phase transformations generated by Eq (25) in a proper way. We will discuss this 
in the following sections. 

Second, let Mc be a differentiable manifold which contains the classical degrees of 
freedom. We shall assume it to be a phase space, and thus it will have an even 
number of degrees of freedom and it shall be endowed with a non degenerate Poisson 
bracket that in Darboux cordinates {R, P) reads 

Third, we let our state space S be the Cartesian product of both manifolds, 

5 := Mc X Mq. 

Such a description has important implications: it is possible to consider each 
subsystem separately in a proper way but it is not possible to entangle the 
subsystems with one another. As long as Ehrenfest dynamics disregards this 
possibility, the choice of the Cartesian product is the most natural one. 

The coordinates that describe the state of the system are: 
The positions and momenta of the nuclei: 

{R,P)eMc. (40) 
We will have 3Nc + 'iNc of these, for Nc the number of nuclei of the system. 
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• The real and imaginary parts of the coordinates of the Hilbert space vectors with 
respect to some basis: 

{q,p)eMQ. (41) 

We will have Nq + Nq of these, for Nq the complex dimension of the Hilbert space 

n. 

5.2. The observables 

Our observables are functions defined on the state space 5* = Mq x Mq. We know from 
our discussion in the case of a purely quantum system that any function of the form 



(29) produces an evolution, via the Poisson bracket, which preserves the norm. In the 



MQCD case, we can easily write the analogue of the vector field (25) by writting: 

Tq = I ® r. (42) 

This is again the infinitesimal generator of phase transformations for the quantum 
subsystem, but written at the level of the global state space Mc x Mq. A reasonable 
property to be asked to the functions chosen to represent our observables is to be 
constant under this transformation. From a mathematical point of view, we shall define 
the set of possible observables, O, as the set of all C°°-functions on the set Mc x Mq 
which are constant under phase changes on the quantum degrees of freedom, i.e., 

0:={fe C^{Mc X Mq)\ TQf = 0}. (43) 

As we shall see later, this choice reflects the fact that, when considered coupled 
together, the nonlinearity of classical mechanics expands also to MQCD. There are 
some important subsets that should be considered: 

• The set of classical functions: these are functions which depend only on the classical 
degrees of freedom. Mathematically, they can be written as those functions / G (9 
such that there exists a function fc G C°°{Mc) such that 

f{R,P,q,p) = fc{R,P). 

We denote this subset as Oc- An example of a function belonging to this set is the 
linear momentum of the nuclei. 

• The set of generalized quantum functions: functions which depend only on the 
quantum degrees of freedom and which are constant under changes in the global 
phase. Mathematically, they can be written as those functions / G (9 such that 
there exists a function /g G C°°(Mq) such that 

f{R,P,q,p) = fQiq,p); T{fQ) = 0. (44) 

We denote these functions as Oq. We have added the adjective "generalized" 
because this set is too large to represent only the set of pure quantum observables. 
These later functions, should be considered, when necessary, as a smaller subset. 



which corresponds to the set of functions defined in Eq. (29). We denote this 
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smaller subset as Oq. An example of a function belonging to Oq is the linear 
momentum of the electrons. 

• A third interesting subset is the set of arbitrary linear combinations of the subsets 
above, i.e., those functions which are written as the sum of a purely classical 
function and a purely quantum one: 

f{R, P, q,p) = fc{R, P) + fQ{m- (45) 

We will denote this set as Oc+q- An element of this set of functions is the total 
linear momentum of the composed system. 

We would like to make a final but very important remark. We have not chosen the 
set of observables as 

V := [fA e C^{Mc X MQ)\fA = (^|A(i?,P)V;)} , (46) 

for A{R, P) a linear operator on the Hilbert space "H depending on the classical degrees 
of freedom because of two reasons: 



It is evident that the set above is a subset of (43) and thus we are not loosing any 



of these observables. But it is a well known property that Ehrenfest dynamics is 



not linear it does not preserve V. We must thus enlarge the set (46). 

We are going to introduce in the next section a Poisson bracket on the space of 
observables. For that bracket to close a Poisson algebra, we need to consider the 



whole set (43). 



It is important to notice that in the set (43) there are operators which are not 
representing linear operators for the quantum part of the system and hence the set of 
properties listed in Section |4] for the pure quantum case are meaningless for them. But 
this is a natural feature of the dynamics we are considering, because of its nonlinear 
nature. 



5.3. Geometry and the Poisson bracket on the classical- quantum world 

As we assume that both the classical and the quantum subsystems are endowed with 
Poisson brackets, we face the same problem we have when combining, from a classical 
mechanics perspective, two classical systems. Therefore it is immediate to conclude that 
the corresponding Poisson structures can be combined as: 

{■,-}:={-,-}c + rH-,-}Q, (47) 

where the term {■,-}c acts on the degrees of freedom of the first manifold and {■,-}q 
acts on the degrees of freedom of the second one. It is a known fact [H] that such a 
superposition of Poisson brackets always produces a well defined Poisson structure in 
the product space. 

The set of pure classical functions Oc and the set of quantum generalized functions 
Oq are closed under the Poisson bracket. The same happens with the quantum functions 
Oq and the set of linear combinations Oq+q- In mathematical terms, what we have is 
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a family of Poisson subalgebras. This property ensures that the description of purely 
classical or purely quantum systems, or even both systems at once but uncoupled to 
each other, can be done within the formalism. 

Once the Poisson bracket on Mc x Mq has been introduced, one can prove, 
in analogy to the full quantum case, that the infinitesimal generator of the phase 
transformations is the Hamiltonian vector field of the identity operator, Tq = {/i, ■}. 
Therefore, a function is a legitimate observable if it commutes with fi. {/n, /} = 0. 



5.4- The definition of the dynamics 

Analogously to the description given in Sections [3] and |4| the dynamics of the mixed 
quantum classical systems can be implemented on 

• The manifold which represents the set of states by defining a vector field whose 
integral curves represent the solutions of the dynamics (Schrodinger picture). 

• The set of functions (please note the differences between the classical and the 
quantum cases) defined on the set of states which represent the set of observables of 
the system. In this case the Poisson bracket of the functions with the Hamiltonian 
of the system defines the corresponding evolution (Heisenberg picture). 

Both approaches are not disconnected, since they can be easily related: 

XH = {fH,-}, (48) 

where we denote by Xh the vector field which represents the dynamics on the phase 
space and by fn the function which corresponds to the Hamiltonian of the complete 
system. 

We can now proceed to our first goal: to provide a Hamiltonian description of 
Ehrenfest dynamics in terms of a Poisson structure. We thus define the following 
Hamiltonian system: 

• A state space corresponding to the Cartesian product Mc x Mq. 

• A set of operators corresponding to the set of functions O defined in Eq. (43). 



The Poisson bracket defined in Eq. (47). 

And finally, the dynamics introduced by the following Hamiltonian function: 

where He is the expression of the electronic Hamiltonian, Mj are the masses of the 
classical subsystem of the nuclei and ipiQjP) is the real-space representation of the 
state ip analogous to \ip) in Eq. (p9|). 



As a result, the dynamics of both subsystems are obtained easily. In the Schrodinger 
picture we obtain: 

k = ^= M-^P, (50) 
dP 
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P = = -gTe.d{{^P{q,p)\H,{R)^{q,p))), (51) 

oR 

qk = h-'^, k = l,...,NQ, (52) 
opk 

Pk = -h-'^, k = l,...,NQ. (53) 
oqk 



This set of equations corresponds exactly with Ehrenfest dynamics. 

Finally, it is important to verify that the dynamics preserves the set of observables 
O, i.e., for any observable /, {/h, /} G O. This can be easily proven since we have 
established that an observable belongs to O if it Poisson-commutes with /j. Thus, as 
fn G O, if we consider an observable / G (9, by the Jacobi identity: 

{/n, {h, /}} = -{/, {/n, fn}} - {fn, {/, fi}} = 0. (54) 



6. A phase space description of the statistics of the Ehrenfest dynamics 

The next step in our work is the definition of a statistical system associated to the 
dynamics we introduced above. The first ingredient for that is the definition of the 
distributions we shall be describing the system with. 

The main conclusion from the previous section is that Ehrenfest dynamics can 
be described as a Hamiltonian system on a Poisson manifold. Therefore, we are in 
a situation similar to a standard classical system. We have seen that the dynamics 



preserves the submanifold Mc x Sq {Sq being the sphere in Eq. (26)), and thus it is 
natural to consider such a manifold as the space of our statistical system. 

We have thus to construct now a statistical system composed of two subsystems. 
We may think in introducing a total distribution factorizing as the product of a classical 
and a quantum distribution. But as the events of the classical and the quantum regimes 
are not independent, the two probabilities must be combined and not defined through 
two factorizing functions. 

Hence we shall consider a density Fqc G O and the canonical volume element dfiqc 
in the quantum-classical phase space, restricted to Mc x Sq: 

d^Qc := d^icd^Q , (55) 

which allows to define the macrosopic average for any observable M E O: 

(M):= / M{R,P,q,p)FQc{R,P,q,p)dfiQc. (56) 

JMcXSq 

In the rest of the section we shall discuss in more detail the ingredients of this definition. 

We start by the volume form. Given the classical phase space Mc with Darboux 
coordinates {R,P) (they always exist locally [H]) we define the volume element in the 
classical phase space by 

diJc:=Y[dPjdRj, (57) 
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which, as it is well known, is invariant under any purely classical Hamiltonian evolution, 
and does not depend on the choice of coordinates. 

We might proceed in the same way for the quantum part of our system. Indeed, as 
we have written the quantum dynamics like a classical Hamiltonian one we could use 
the canonical invariant volume element describe in the previous paragraph d^q on Mq. 
However, new complications appear in this case as in the quantum part of the phase 
space we restrict the integration to Sq and then we look for an invariant volume form 
in the unit sphere, rather than in the full Hilbert space. 

The new ingredient that makes the restriction possible is the fact that all our 
observables and of course, also the Hamiltonian, are killed by F. This implies, as 
discussed before, that fi is a constant of motion. Using this we decompose 

d^Q = dfidQq, (58) 

which is possible on \ {0}. The decomposition is not unique, but the restriction 
(pullback) of dQg to the unit sphere gives a uniquely defined, invariant volume form 
on Sq that we will denote by dQq. The decomposition is nothing but the factorization 
into the radial part and the solid angle volume element. 

Example 2 For the simple case we studied above, where T-i = and hence Mq = 
and Sq is a three dimensional sphere, the volume element above would be: 

dVlq = qidq2dpidp2 - q2dqidpidp2 + Pidqidq2dp2 - P2dqidq2dpi. (59) 

Now, we finally put together the two ingredients to obtain the invariant volume 
element dfj,qc in eq. (55). 



We discuss now the properties we must require from Fqc in order for Eq. (56) to 
correcly define the statistical mechanics for the Ehrenfest dynamics. 

• The expected value of the constant observable should be that constant, which 
implies that the integral on the whole set of states is equal to one: 

Fqcdi^qc = 1- (60) 

McXSq 



• The average, for any observable /a of the form (46) associated to a positive definite 
hermitian operator A should be positive. This implies the usual requirements of 
positive probability density in purely classical statistical mechanics. 

It is simple to go from this phase space description to the formulation of a density 
operator p for the quantum system and make the construction closer to the usual 
description of Quantum Statistical Mechanics. The density operator is defined as follows: 

p{R,P):= [ dpq{q,p)FQc{R,P,q,m{q,P)){mP)l (61) 

where \4'{q,p)){ip{q,p)\ is the projector on the quantum state parametrized by the pair 
{q,p). The average value of a pure quantum magnitude A, represented by the quadratic 
funtion f^ can be computed as: 

{A){R,P):=Tt{p{R,P)A)= I dpq{q,p)FQciR,P,q,p)fA{q,p). (62) 

JSq 
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To obtain the corresponding total average, it is required to integrate the resulting 
object on the classical phase space: 

(A) := I dfic{R,P){A){R,P). (63) 
Jmc 

So far we have discussed general considerations about the statistical mechanics 
of our system. Now we should include the dynamics in the game. Given that we have 
formulated the Ehrenfest dynamics as a Hamiltonian system, it is well known [37] which 
is the evolution of the density function Fqq- It is given by the Liouville equation 
d F 

-^ + {fH,FQc} = 0, (64) 

where, in the derivation of the equation, it is an essential requirement the invariance of 
the volume form under the evolution of the system. 

The equilibrium statistical mechanics is obtained by requiring Fqc = 0, which 



making use of (64) is equivalent to 

{fH,FQc} = 0, 

or, in other words, Fqc should be a constant of motion. An obvious non trivial constant 
of motion is given by any function of fn- The question is if there are more. To answer 
this we should examine closely our particular dynamics. In several occasions we have 
mentioned the non linear character of a generic Ehrenfest system. This is an essential 
ingredient, because due to this fact we do not expect any other constant of motion that 
a function of the Hamiltonian itself. This is in contrast with the case of linear equations 
of motion, where the system is necessarily integrable. 

The question of the existence of more constants of motion and, therefore, more 
candidates for the equihbrium density Fqc can be elucidated with a simple example 
that we discuss below. In it we show the presence of ergodic regions (open regions 
in the leave of constant energy densely covered by a single orbit). This rules out the 
existence of constants of motion different from the constant function in the region (i.e., 
functions of fn)- 

Example 3 In the following example we will study a simple toy model in which the 
coupling between classical and quantum degrees of freedom gives rise to chaotic behavior 
and the appearence of ergodic regions. 

The system consists of a complex two dimensional Hilbert space Mq = and a 
classical 2-D phase space where we define a 1-D harmonic oscillator. Using coordinates 
{le, 6) for the classical variables (action-angle coordinates for the oscillator) and G 
we define the following Hamiltonian 

fH = h + + ecos(e)a^|^), 

with 0"^, the Pauli sigma matrices. 
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Figure 1. The plot shows a single orbit of the Poincare map of Eq. (661 at 6* = 0. 
The angle cj) G '"'] is plotted in the horizontal axis and G [0, 1] in the vertical 
one. We have taken e = 0.8. The energy determines 1$ and its actual value does not 
affect the dynamics of (h and Ij,. 



We parametrize the normalized quantum state by 

h. 



1^) = e*" 



(65) 



One can now easily check that and are canonical conjugate variables. In these 
variables the Hamiltonian reads 



fH = h + ll + el^^l - II cos(^) cos(0), (66) 

where e measures the coupling of the classical and quantum systems. 

In the limit of vanishing e the system is integrable and actually linear in these 
coordinates. However for non vanishing e the model becomes non linear and, as we will 
show below, regions of chaotic motion emerge. 

In order to understand the behaviour of our system it is useful to study its Poincare 
map (see Ref. for the definition). To this end we take the transversal (or Poincare) 
section atO = and, taking into account the conservation of energy, only two coordinates 
are needed to describe the map, we have chosen the quantum variables I^ and 

By examining the plot of Fig. [IJ one sees that the orbit densely fills a large region 
of the Poincare section and therefore any constant of motion should be the constant 
function in that region. In the presence of such an ergodic evolution the only constants 
of motion are functions of the Hamiltonian. 

As a consequence we can claim that, generically within the Ehrenfest dynamics, 



the only functions which commute with the Hamiltonian in (64) are those which are 



function of the Hamiltonian fu itself. As the Poisson bracket is skew-symmetric, the 

I In our construction of the Poincare map we represent the value of (j) and whenever the periodic 
coordinate 9 takes the value 6 = Q. The successive values of (f) and I^, for some initial condition, 
form our orbit. As stressed in [4r the Poincare map is a very useful tool in the context of non linear 
dynamics as many properties of the full system like ergodicity, chaos or regular behaviour, may be 
inferred from the same property of the Poincare map. 
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property follows trivially. Ergodicity helps us to assert, roughly speaking, that in general 
no other function in O will be a constant of the motion. 

Thus, at this stage, we can finish our construction with the equilibrium distribution 
associated to Ehrenfest dynamics. Taking into account the ergodicity of the dynamics, 
we can impose the equal-probabilities condition to the configurations of an isolated 
Ehrenfest system, which leads to the microcanonical ensemble, (see Ref. ^3]). 

In the case in which this system is weakly coupled to a bath, if we assume that 
the dynamics of the system plus bath is ergodic as well, such that system plus bath are 
microcanonically distributed, then the statistics of the Ehrenfest system are given by 
the well known canonical probability density (see, for example, Refs. [3HI SI]): 

F^aiR, P, q,p) = Z-'e-^f", Z = [ e-^^^df^Qc , (67) 

where proportional to the temperature of the subsystem, is the parameter governing 
the equilibrium between the system and the bath. An information-theoretic approach to 



the equiprobability in the microcanonical ensemble and to Eq. (67) may be seen in Ref. 
[15] or in Refs. |l6l HTj for a much more modern and mathematically sound exposition. 
We also strongly recommend the recent work of P. Reimann [IHl HQ] for a foundation of 
Statistical Mechanics under experimentally realistic conditions. 

7. Application: Nose formalism for Ehrenfest systems 

We now proceed to extend the non-stochastic method proposed by Nose pTl [TS] to our 



Hamiltonian (49). Consider a system defined on Mc x Mq as before and an extension 
of that system including two new classical degrees of freedom {s,Ps), what produces a 
total state space 

Mcx MqX 

In this case, in order to include the quantum degrees of freedom, we must take into 
account the problem of the deformation of the domain (see the change of coordinates 



of Eq. (77)) of integration we dealt with in the previous section. Thus we must rewrite 



the Hamiltonian given in Eq. (49) in order to extend the domain of the quantum part 



from the sphere to the whole Mq. In order to do this, we consider: 

=:H'^iP) + eH^iq,p,R), (68) 

Now we transform this Hamiltonian into the "extended" one, by introducing a 
dependence with respect to the new degrees of freedom: 

2 

H{R,P,q,p,s,p,) := H^iP/s) + eH^{q,p/s,R) + ^ + gkT\og{s) . (69) 

This is very similar to the method proposed by Nose [T71 [18] for thermostating 
classical systems. The Nose technique is basically an algorithm to produce the 
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(74) 



correct equilibrium distribution and, therefore, the interpretation of the coordinate s as 
describing the behaviour of some 'thermostat' is only metaphorical. In the same sense, 
the particular form of the last two terms in the above Hamiltonian is just the one needed 
for the theorem that we mention in the following lines to work. 
The extended equations of motion read now: 

R = P/{Ms^), (70) 
P = -grad^ (^eH-{q,p/s,R)^ , (71) 

<ik = (^1^) i^Pk^iQi + ■■■ + Hp^pnqPNq/s- eH-{q,p/s,R)pk), (72) 

Pk = -^ji^{Hq^q^qi + . . . + Hg^p^^pNg/s- eH-{q,p/s,R)qk), (73) 

. Ps 

p,^s-^J2f + Y.^'-'''sI<T, (75) 

k j 

where k = 1, . . . , Nq. 

Now we generalize the standard result of Nose [T71 dS] , by proving the equivalence of 
the microcanonical distribution of the extended system and the canonical distribution 
of the original one. The equations (70 75) are the necessary ones to implement the 
microcanonical scheme once we assume the ergodic hypothesis for the extended system. 
To prove the equivalence with the canonical distribution of the original system, we take 
the Hamiltonian above, and we consider the corresponding microcanonical distribution: 

Z= dficdHQdpsdsS{H{R,P,q,p,s,Ps) - E), (76) 

J McxMqxR'2 

where ducd^Q is the volume element on Mc x Mq. 

If we consider the following noncanonical transformation on Mc x Mq x : 

I^'=R\ P;=^, q" = q\ p[ = ^, s' = s, p:=Ps, (77) 
s s 

the partition function reads (we use the same symbol again for the transformed variables 

but include the jacobian which is equal to s^): 

Z = / dficdfiQdpsds 

J McxMqxR^ 

6 (^if^(P) + CH^q.P, ^) + ^ + gkThgis) - E^ , (78) 

where D is the number of degrees of freedom of the system (the sum of classical and 
quantum degrees of freedom). 

As the original integral domain includes the complete Mq, and therefore a rescaling 
of momenta does not modify it, we can use the equivalence 



5{f{s)) = 6is-so)/\f'iso)l 
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where 



and write: 



So = exp 



/7^(P) + e^,e(g,p,i?) + g-E 



gkT 

dncdfJ,Qdpsds 



s — exp 

Integrating in the variable s, we get: 



H^(P) + eH^{qJ,R) + l^-E 



gkT 



(79) 



gkT 

j 



dps exp 



D + 1 



gkT 
dficdfiQ exp 



E 



Pt_ 
2Q 



McxMq 



D + 1 



iH''{P) + eH4q,p,R) 



gkT 



(80) 



Then if we choose g = D + 1 the exponential takes the desired form. Besides, we can 
notice that when calculating the average of any function M G (9 (where O corresponds 
to Eq. (43) and therefore M does not depend on the extended variables), the quantity: 



1 



{D + l)kT 



dps exp 



1 
kf 



E 



P\_ 
2Q 



factorizes. If we also assume that the function M does not depend either on the norm of 
the quantum state (as it is the case for the quantum observables) and we decompose the 
volume element dfiq as Eq. (58), the radial part of the quantum integral also factorizes 
and cancels the analogous term arising from the partition function. 
Thus we can write 



(M) = 1/Z 
= 1/Z 



dficdfidnQdpsdsM{R, P, q,p)diH{R, P, q,p, s,ps) - E) 

1 

D + l)kT 



dps exp 



1 

kf 



E 



2Q 



d/ix 



dficd^QM{R, P, q, p) exp 

MqXSq 



1 

kf 



H^{P) + eH^{q,p,R) 



where 



Z 



1 



dps exp 



{D + l)kT 

dficdi^Q exp 



1 

kf 



E 



2Q 



dfi X 



1 

kT 



H''iP) + eH^iq,p,R) 



'McxSq 

Hence we have proved that the microcanonical distribution of the extended dynamical 
system created by Eq.(69) is equivalent to the canonical distribution created by the 
Hamiltonian 

H{R,P,q,p} = H''{P) + eH4q,P,R)- 
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8. Conclusions and future work 

In this paper we liave constructed a rigorous Hamiltonian description of the Ehrenfest 
dynamics of an isolated system by combining the Poisson brackets formulation of 
classical mechanics with the geometric formulation of quantum mechanics. We have 
also constructed the corresponding statistical description and obtained the associated 
Liouville equation. Finally, after verifying numerically that the Ehrenfest dynamics is 
ergodic, we justify the equilibrium distribution produced by it. 

It is important to keep in mind the restriction to finite-dimensional quantum 
systems that we fixed as starting point, since our first goal is to apply our construction to 
Molecular Dynamics. From a theoretical point of view, though, it is very interesting to 
perform the same type of analysis for infinite dimensional quantum systems (note that 
the Ehrenfest dynamics as formulated in [U [5] consists of classical nuclei and an infinite- 
dimensional Hilbert space for the electronic subsystem). Formally, the construction 
can be defined in the same way, once we adapt the language to infinite dimensional 
manifolds. But the nonlinearity of the resulting dynamics imposes very challenging 
problems at many levels, for instance at the very existence of the solutions of the 
Hamiltonian dynamical fiow. Another interesting problem would be to recover, within 
the geometrical formalism, the classical-quantum system as a suitable classical limit for 
the nuclear degrees of freedom of an original completely quantum system (similarly to 
what is done in [U [5] in the standard approach). We hope to address these problems in 
the future. 



The definition of the Hamiltonian in Eq.(49) and its associated canonical 



equilibrium distribution Eq. (67) allows us now to use (classical) Monte Carlo methods 



for computing canonical equilibrium averages in our system, given by the expression Eq. 



(56). But this formalism also allows to easily introduce a Nose extended system scheme 
to compute averages in the canonical ensemble by performing constant energy Molecular 
Dynamics simulations. This is a very important feature of the final setting, because the 
situations where the Ehrenfest approach is necessary (i.e. those where the electronic 
energy gap is small and the Born-Oppenheimer approximation is not sufficient) the 
effect of the temperature can be relevant. A method for the simulation of the canonical 
ensemble is thus of the greatest importance. 

Alternatively, if we wish to perform MD simulations using stochastic methods, we 



should first construct the Langevin dynamics associated with our Hamiltonian Eq. (49) 



and equations Eqs. (50)-(53). To develop this program, we also need an analogue of the 



Fokker-Planck equation associated with those Langevin equations and then we have to 



check whether its solution at infinite time approaches Eq.(67). Presently we are working 



on this point [50], i.e., on the extension of our formalism to an associated stochastic MD. 
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